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ABSTRACT 

We have examined 93 SSC's during the four year period from 1968 to 
1971 at and near the peak of the solar activity cycle. Of the 93 SSC's 
81 could be associated with solar activity, such as solar flares and 
radio bursts of Type II and Type IV. The mean propagation speeds of 
these flare-associated events ranged from 400 to 1000 km/sec with an 
average speed of 600-700 km/sec. Disturbances associated with 48 of the 
SSC's have been studied in detail using the corresponding interplanetary 
(IP) magnetic field, and plasma data when they were available. We found 
that 41 of the 48 disturbances corresponded to IP shock waves, and the 
remaining seven events were tangential discontinuities. Thirty percent 
of the IP shocks had thick structure (i.e. the magnetic field jump across 
the shock occurred over a distance much greater than 50 proton Larmor 
radii). Also given is a statistical study of the gross geometry of a 
"typical" or average shock surface based on multiple spacecraft sightings 
ana their relative orientation with respect to the solar flare and/or 
solar activity. By considering the orientations of 22 well-determined 
shock normals in relation to the positions of the parent flares on the 
solar disk, which is a modification of a method given by Taylor (1969), 
it is suggested that a typical shock front propagating out from the sun 
at 1 AU has a radius of curvature on the order of 1 AU. Also given are 
some general properties of oblique IP flare-shocks. 


Introduction 


It is generally believed that most IP shock waves originate at or 
near the sun, in particular from a solar active region (Gold, 1955; 
Hirshberg, 1968; Hirshberg et al., 1970; Hundhausen, 1970; and Hundhausen 
et al., 1970). Some theoretical models of the propagation of flare- 
associated shock waves in an ambient solar wind based on similarity 
solutions of the hydrodynamic equatixsra, were developed by Parker (1961), 
Simon and Axford (1966), Dryer (1970), and Korbeinikov (1969). Hundhausen 
and Gentry (1969a, b) used numerical solutions of hydrodynamic equations 
to simulate the propagation of flare-associated disturbances. De Young 
and Hundhausen (1971) found that, even for a blast confined to a cone of 
half-angle equal to or less than 15°, the shock front upon reaching 1 AU 
becomes quasi-spherical centered at approximately 0.5 AU. For a model 
of flare-associated IP shocks, shocks are expected to form in the 
vicinity of the sun and propagate outward with a shock thickness on the 
order of a few proton Larmor radii during most of their passage through 
IP space. Spherical symmetry of the shock surface may be disturbed due 
to IP inhomogeneities (Heinemann and Siscoe, 1973) or abrupt discontin- 
uities (Lepping and Burlaga, 1973). 

Hirshberg (1968) derived an average shape of an IP shock front 
from a statistical study of the magnitudes of SSC's. She concluded 
that flare ejected plasma is emitted on a broad front but with consid- 
erable departure from heliocentric spherical symmetry. Taylor (1969), 
utilizing IMP 3 observations, found that the shock front at 1 AU can 
be well approximated by a spherical surface whose curvature is ^ 0.75 
AU and centered at ~ 0.5 AU. 


In this paper we present an observational study of 38 flare- 
associated shocks. The shock speeds and normals have been computed 
accurately for 22 of these cases through the use of multiple spacecraft 
(S/C) observations; the speeds and normals for the remaining 16 shocks 
are also obtained. 

The Experiments 

IP data have been collected from eight S/C (Explorers 33, 34, 35, 

41, 43, Pioneers 7 and 8, and 0G0 5), whicfy taken together cover the four 
year period of interest (1968 through 1971). We are primarily dealing with 
IP magnetic field data and the positions of pertinent S/C during the time 
of onset of the events. Only in a few cases were the detailed IP plasma 
data available to us; we used them mainly to check results obtained from 
the multiple S/C method and to obtain local plasma bulk speeds. Table 1 
shows the S/C and the associated principal investigators for the experiment 
from which we obtained this data. Many of the events have been observed by 
more than one S/C. However, in some intervals only the magnetic field 
data were available. Hence, only a limited study could be made for those 
periods. 

The Method of Study 

Sudden commencements (SSC's) identified by 20 or more magnetic 
observatories were selected for the years 1968 to 1971 inclusive. Then 
for a given event the IP magnetic field data from the experiments on the 
S/C listed in Table 1 were examined. For the events for which the IP 
data were available the associated discontinuities In the magnetic field 
were selected. When the events had been observed by more than one S/C, 




Table 1 


Interplanetary Data Source 


Principal Investigators 


Magnetic Field Experiments 


Plasma Experiments 


N. F. Ness, GSFC 


H. Bridge, M. I. T. 


K. W. Ogilvie, GSFC 


S. J. Bame, Los Alamos 
and K. W. Ogilvie, GSFC 


Explorer 35 


N. F. Ness and F. Marlani, Rome 
P. J. Coleman, UCLA 
C. P. Sonett, Ames 
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multiple S/C methods could be used to find the normal and speed of these 

discontinuities. In orde. to find an accurate normal and speed for a 

given event, we used the following procedure: 

First, the magnetic field data were selected for each event. The 

average values and their standard deviations were computed for both sides 

of each discontinuity (here discontinuity means shock or tangential 

discontinuity). The time intervals over which the averages were taken 

depended on the behavior of the fluctuations of the vector quantity B 
in the vicinity of the discontinuity. Time intervals were chosen over 

which this quantity displayed relatively steady behavior; these were in 

general 3 to 5 min long. Then, the shock normal n g was computed for each 

S/C observation using the magnetic coplanarity theorem (Colburn and Sonett, 

1966). The uncertainty of n g is closely related to the fluctuations of 

the magnetic field in the analysis interval. However, for most IP events 

the computed n g 's have an uncertainty within Z 20°. If the discontinuity 

was a tangential one, the normal to the discontinuity plane should be in 

Bj x ^ 

the direction n where n„ ■ — . The direction n,. was also 

' ' | (B X B,)| ' 

computed for each observation, where and are the average magnetic 
field before and after the discontinuity, respectively. The comments 
above concerning the expected error in n g hold for n t ; however, in 

A ^ -a 

general the error in n^ should be smaller, since (B^ - B^) is not involved. 
Then the normals n g and n t were computed in every case regardless of 
whether the actual event corresponded to a shock or tangential discon- 
tinuity, and ideally they should be at 90° with repect to each other. 

Since in practice n g and n^ can be distinguished from each other with 


. 


/ 
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only an error of £ 30° (=/? 20°), the ideal 90°, being a factor 3 
greater, generally enables highly reliable differentiation. Choosing 
which of the two normals was correct is another matter and is discussed 
below. 

Secondly, multiple S/C methods, which were also used, are now 
described. If there were two S/C observations available, the relative 
position vector of these two S/C, AR - (R 2 “ R^X is related to the 
propagation speed of the discontinuity, V^, and to the time difference, 
At, between observations of the discontinuity at these two S/C in the 
following way: 

AR * n d - V d At (1) 

where and R^ are the position vectors of the two S/C and n d and 
are the normal and the local speed of the discontinuity, respectively. 
Note that we do not yet specify the type of discontinuity, i.e., shock 
ot tangential. This calculation is based on a geometric configuration 

only— and kinematic assumptions to be discussed — but is independent 

A 

of specific identification. Also since the normal component of the 
magnetic field is continuous across any type of discontinuity, then 
AB (^(B^ - B^)) is parallel to the surface of the discontinuity, i.e., 
the condition 

AB * n . - 0 (2) 

a 

holds for any type of discontinuity. It is important to note that on a 
scale of the order of 0.01 AU a shock front usually remains planar and 
propagates at a constant speed (Ogilvie and Burlaga, 1969; Chao, 1970), 
and similarly for tangential discontinuities on the order of 0.002 AU 


r i * . //i f 




(Burlaga and Ness, 1969). In all of the computations done for this 
study Involving more than one S/C, we have assumed In fact that the 
discontinuities can be approximated as planar and propagate at a constant 
speed within the distance of the S/C separations (% 0.004 AU) . However, 

-4 

the magnetic field can change over all scale lengths. That Is, the AB's 
are not necessarily equal to each other for the two S/C observations, 
even though the normals remain the same. Hence, It was often useful to 
provide two equations of the type Eq(2), one each for the two S/C 
observations. Then solving Eqs(l) and (2) together yielded 9nd n^ 
for the discontinuity. 

If three S/C observations were available, two equations of the type 

Eq (1) and any one equation of the type Eq(2) were enough to determine V d 

and n,. In the case of four S/C observations having been available, three 
d 

equations of the type Eq(l) were used to determine n^ and V^. This case 
was completely kinematically determined from the geometric configuration 
and onset times alone. 

For a few cases that were well known to be shocks, and when plasma 
data were available, one of the best-fit methods using a subset of the 
shock conservation equations was emp >yed to obtain n g (Chao, 1970; 

Lepping and Argentiero, 1971). 

The normal n., computed from the multiple S/C method, was compared 
d 

with both r» t and n g , obtained from the single S/C magnetic field measure- 
ments. And when available the solar wind speed was compared to the 
estimated discontinuity speed divided by the radial component of the 
unit normal. It is, in general, possible to differentiate a shock wave 
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from a tangential discontinuity by these comparisons. Well known shock 
and tangential discontinuity signatures in the magnetic field data 
(Burlaga, 1968) also may be used as a guide in the discrimination. 

By combining the single and multiple S/C methods, it was possible in 
general to lessen the uncertainty of the estimate of the normal of the 
discontinuity, yielding final errors usually of approximately + 10°. 

The preceeding scheme and assumptions form the basis for the method used 
to obtain the identification and kinematic properties of th discontinuity. 
The Observations 

For the active period of the solar cycle from 1968 to 1971 inclusive, 
we selected 93 world-wide SSC's which have been reported by more than 20 
geomagnetic stations each (Solar-Geophysical data, published by ESSA, U.S. 
Department of Commerce). Then the solar activity which occurred one to four 
days before the SSC's were examined. The SSC's usually could be associated 
with a flare of importance IB or stronger and radio bursts of Type II and 
Type IV. On an average the solar activity association of these events can 
be made with a reasonable degree of certainty, as we will demonstrate in 
a later section. 

Of the 93 SSC's 81 could be associated with solar activity. Eighty- 
five percent of these associations included radio bursts of Type II or 
Type IV. Therefore, we believe our statistics of these 81 events are 
significant. 

The IP magnetic field data from the S/C listed in Table 1 were 
available for 48 of the 93 SSC events. From these IP data our analysis 
shows 7 of the 48 events were tangential discontinuities and the remaining 


41 were shock waves. However, only 38 of these 41 shocks had adequate 
IP data to obtain shock normals and speeds. 

Chao (1973) has shown two examples of IP shock-like discontim ities 

which had a "shock" transition region of a thousand proton Larmor radii 

or wider. These two events steepened into shock waves in the vicinity of 

the earth. These two "shocks" can be associated with solar flares and 

with radio bursts of Type il and Type IV. We examined the thickness 

of the magnetic field transitions of the 38 shocks and found that 307. of 

the "shocks" had a transition region larger than 50 proton Larmor radii, 

Rp. In most of the thick-structure events, the transition regions were 

more than a few hundred R in width. We do not claim that the events with 

P 

a thick transition zone are fully formed shocks. We would like to suggest 
that what we arc observing in these ca«es is the formation of shocks. 

Out of the 38 shocks, 22 shock normals and speeds were computed 
accurately, and their solar activity associations were relatively 
reliable. We will study IP shock correlation with solar activity using 
these 22 events. Table 2 gives a summary of the results discussed in 
this section. 

Results 

If we assume that the 81 flare-associated events represent distur- 
bances tuch as shock waves propagating over 1 AU from the flare site to 
the earth's vicinity, then the average transit velocity of such a distur- 
bance can be computed using the difference between the occurrence of the 
flare and the onset of the SSC at earth. A histogram of these mean speeds 
is given in Figure 1. There is a peak at 600 - 700 km/sec. The 
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Table 2 


Some Statistics of SSOAssocltted Events 

Period of Study: January 1968 to December 1971 

93 world-wide SSC's were selected, 

81 of the SSC's were associated with flares and in some cases radio 
bursts (Type II and/cr IV). 

48 hnd available interplanetary (IP) magnetic field and in some 
cases plasma data for the associated IP events. 


Analysis of 48 IP events yielded: 

7 tangential discontinuities 

f 3 had unknown normals and speeds. 
38 had known normals and speeds: 


41 shocks ^ 

l 


26 were thin structures 
12 were thick structures 


22 of the 38 shocks had accurate shock normals and speeds estimated 
and all 22 were flare associated. 


distribution in the upper and lower speed portions may not be reliable 

due to the difficulty of selection for the extremes. However, the peak 

of the distribution is significant. This mean speed of *00-700 km/sec 

implies that the shock speed is not on an average very much greater than 

the solar wind speed, i.e., the shock is of intermediate strength. 

The shock normals of the 38 shocks are shown in Figure 2. The normals 

are given in solar ecliptic coordinates where 0^ is the latitude angle 

and <f> is the azimuth angle measured in the ecliptic plane. The dotted 

arrows represent those events with larger uncertainties than those with 

solid arrows; the latter set (22 shocks) usually correspond to multiple 

spacecraft observations and are well determined. These normals are rather 

symmetrically distributed with respect to the sun-earth line. 

For those 15 of the 38 shocks for which plasma data were available 

the dependence of the shock speed in the solar wind frame of reference 

on the local solar wind speed in the pre-shock state was studied. Figure 

3 shows the local shock speed in the solar wind frame, W, versus the 

solar wind speed V . The figure shows that the pre-shock solar wind 
ow 

speeds occur in the range 350 to 420 km/sec, which corresponds to the 
most probable solar wind conditions. W, however, ranges from 40 to 300 
km/ sec. 

Next, we compare the mean shock speeu with the local shock speed 

as computed by the multiple spacecraft method. Figure 4 shows a plot 

of the mean transit velocity <Vg> versus the local shock speed V g . 

Notice that most of the events lie below the diagonal line. That la, the 

mean speed <V > is larger than V for most events. This shows that, in 
S b 

general, flare-produced shocks slow down during propagation over 1 AU. 
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The orientation of shock normals and their positions relative to the 

"arent flares can be studied to obtain an average shock front. Individual 

locally-determined shock surface orientations may severely deviate from 

this front (Lepping and Chao, 1972). The 22 shock events (for which 

normals have been accurately computed), their normals, and information 

about the identification of the probable parent flares are listed in the 

top section of Table 3. The remaining 16 of the 38 shocks are listed at 

the bottom of the table. The first two columns give the dates and times 

of the SSC events, respectively. The third column gives the Cod» number which 

identifies each event. The fourth through eighth columns give the following 

for each flare: the time of onset, position in latitude, importance, and 

the time of ontet of radio bursts of Type II and Type IV, respectively. The 

next three columns give the shock normal in solar ecliptic coordinates 0 g 

and <t> (<t> = 0° in the solar direction), and the local speed computed from 

s s 

the multiple S/C method, respectively. The average sun-earth transit 
speed is given in the last column. 

These 22 events are related to their parent flares in Figure 5 
according to a method by Taylor (1969), and recently repeated by Bavassano 
et al. (1973) using the Pioneer 8 data. Figure 5 is a plot of the 
orientations of the 22 shock surfaces in the ecliptic plane at the appro- 
priate heliocentric longitude relative to the flare. The longitudes are 
specified as seen from the earth. For example, the flare associated with 
the SSC of November 20, 1968 occurred at 90°W on the solar disk and thus 
is plotted at 90°E of the flare as shown in the figure (Code no. 19). 

This figure shows that the average shock surface in the ecliptic plane 
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near the earth's orbit lies on a circle centered at the sun with a radius 


of 1 AU. A few events have been observed near the west limb but these 
also lie well on the circle. The mean deviation of these shock surfaces 
from the circle is 5° (+ 9°), i.e. << *>> = 175° for the average normal, 
with a standard deviation (a) from the mean surface of about 22°. The 
mean deviation of the remaining 16 shock surfaces from the circle is 3° 

(+ 18°) and a % 36°. This is somewhat different from the results given 
by Hirshberg (1968), Taylor (1969) and Bavassano et al. (1973) which 
suggest smallar radii of curvature by approximately a factor of 2. 

Figure 6 shows the orientation of the 22 shock surfaces in the 
meridional plane at the appropriate heliocentric latitude relative to the 
flare position. The shock surfaces are clumped above and below the 
ecliptic plane at about + 20°. The figure shows that there is a tendency 
for the average shock surface in the meridional plane near the earth's 
orbit also to lie on a circle centered at the sun with a radius of 1 AU. 
The mean deviation of these shock surfaces from the circle is 0.6° (+ 10°) 
with a a of about 24°. For the remaining 16 shocks the corresponding 
values are 3° (+ 20°) and 39°, respectively. This implies that the 
average shock surface does not deviate much from a spherical shape in 
the meridional plane. 

Using the 22 flare-associated shock waves, we obtained statistics 
about the correlation between shock waves, flares and radio bursts of 
Type II and Type IV. Table 4 gives the probability of solar activity 
being associated with interplanetary shock waves. We have divided the 
solar activity into three classes, namely: flares, radio bursts of 
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Table 5 


Flare and Radio Bursts-Associated Interplanetary "Random" Shock Waves 


Type of Associations Percentage of IP Random Shocks 


Flare, Type 

II and IV 

207. 

7 lare, Type 

II 

407. 

Flare, Type 

IV 

307. 

Flare only 


507. 

Flare, Type 

II or IV 

507. 
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Type II and Type IV, which are believed to be most strongly related to 
interplanetary shock waves (Wild, 1970). From this table one can see 
that, of radio bursts, Type II have the stronger correlation with IP 
shocks. Note that in general one can always find a flare-association 
with a given IP shock, because theiC are so many flares reported over 
any four day period. On the other hand, if the radio data are included, 
it is seldom possible to make more than one association with a given IP 
shock. This table shows that we apparently do have a meaningful associ- 
ation with IP shocks and solar activity. 

In order to test our method of associating solar flares with the 
observed IP shocks, we artificially generated random IP shock times and 
repeated the association procedure as if the times were real shock onsets 
(or SSC's) at 1 AU. The day of occurrence was generated by computer 
assuming a uniform distribution and an expected rate of 2 shocks per 
month for 2 years (1968-1969). The hour of day was derived by throwing a 
die so that resulting "onset times" had a quantization of 4 hours dura- 
tion. (Over a large number of throws the value showing on the die is 
expected on average to be 3-1/2: in our case of 48 throws it was 3.46.) 
Using these times we attempted to associate solar flares, Type II-, and/or 
Type IV-solar bursts just as we had done in the cases of the real IP 
shocks. For all but two random shocks we were able to find some flare 
association. For these 46 cases a histogram of mean speeds was produced 
as shown as a solid curve in Figure 7, corresponding to Figure 1 for the 
real cases. In some instances radio bursts were associated, and the 
percentage of these associations is shown in Table 5, corresponding to 
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Table 4 for the real cases. For some of the associations it was 
necessary to ignore the more impressive solar flare, by importance 
designation, in favor of a lesser flare because of the existence of radio 
bursts, Type II and/or IV at the time of the lesser flare. This was also 
done in the real cases but not as frequently. Firstly, we point out that 
the histogram in Figure 7 (solid line) shows a skewed distribution with 
respect to the histogram in Figure 1, which is shown in Figure 7 as a 
dashed curve, properly scaled, for comparison. Secondly, the most probable 
speed is shifted to an unrealistically low value of ~ 500 km/sec, instead 
of 650 km/sec for the real cases. (The average post shock solar wind 
speed alone over 1 AU is expected to be greater than 500 km/sec). Thirdly, 
and most important, as Tables 4 and 5 show, the percentage of occurrences 
of solar radio bursts in any combination is much lower than in the real 
cases. In particular compare Flare, Type II or Type IV for random (50%) 
to real cases (85%). Also the flare importance designation was necessarily 
lower, on an average, for the random case associatlor*. Lastly, we point out 
the difference of the mean speed averages but the similarity of the ct's, 
as Table 6 shows. We conclude that the method of using an IP shock-solar 
association based on Flares, and Type II and IV bursts when present, for 
this statistical study yielded reasonably trustworthy results. But we 
caution that ■<ne must be exceedingly careful about individual shock 
studies based on such a method. 

Some Flare-Shock Properties 

In this section we list and discuss those properties of oblique IP 
flare-shocks which depend explicitly on magnetic field quantities and 
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Table 6 


Characteristics of Mean Speed Distributions 



No. of Events 
(N) 

Ave + Error 

a 



in km/sac 


Random Test 

46 

660 + 70 

240 

Real Cases 

82 

790 + 50 

230 

Difference 


130 + 90 

10 
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shock normals only, based on the shocks in Table 3, top. In particular 
we are concerned with giving a general characterization of typical 
magnetic field jumps across such shocks. And we are interested in the 
expected importance of thermal anisotropy on estimating shock parameters 
and normals using a fitting scheme based on shock conservation equations 
(Lepping and Argentiero, 1971). Also we wish to check the accuracy of 
the estimated normals presented in this paper. Only the most accurate 
set, according to the shock normal estimates, of the 22 shocks will be 
examined. These shocks will be characterised by the quantities defined 
in Table 7. In the table column- number referes to the respective column 
in Table 8, where the results are presented. The parameter R in column 
7 provides a means of measuring the importance of anisotropy (Lepping, 
1972). That is, for 0 £ R £ 1/2 and for expected anisotropy condi- 
tions in the vicinity of the interplanetary shock (Hundhausen et al., 1967; 
Chao and Goldstein, 1971) departure from the assumption of thermal isotropy 
is in genral unnecessary (Lepping, 1972). Below we show that indeed R 

— ♦ A 

£ 1/2. Concerning column 9, AB • n = 0 must be satisfied across the shock; 
checking it gives the degree of consistency of estimated normals with 
this constraint, where 6=0° indicates perfect agreement. S and S q ere 
arbitrarily defined measures of the "probable strength" of the shock when 
plasma data is unavailable; ^2^1’ * s a s ^ ra ^^ ar indicator. All other 
quantities in the Table 7 are self-explanatory. 

These quantities were calculated for the 22 shocks and the results 
are presented in Table 8. In obtaining these results the following 
rules were applied. Only data from one S/C were used, a ^ 10°, 6 *£ 10®, 


Table 7 


Definitions of Characteristic Quantities 
of 22 Flare-Shocks 


Column 


1 

2 

3 

4 

5 

6 

7 

8 


9 


10 


11 


Quantity 

Code number refers to an event in Table 3. 


S/C = spacecraft's data that was used. 


l * 4- (n, "Y 
9 = * (n, B ) J 


(B 1§ B 2 ) 


B. = 
8 

a 


where n is the shock normal 
and 1 and 2 are pre- and 
post-shock, respectively. 


8^ + a 

R = tan fl^/tan 

F-/F. where F. = V B ? + 
2 1 i ix 


~2 T~ 

K + B j 
iy is 


1,2 


f 


AB 


B 2 - B, 


6 ■ 90° - cos" 1 (n *|^||). 


\ f 2 + V 


for F^ see 8 above. 


F - F 
2 1 


for see 8 above. 


Table 8 


.Characteristics of 18 Flare-Shork* 


Code 

S/C 


B 2 

a 

+ a 

R 

V F i 

6 

S 

S 

1 

33 

61 ° 

85 ° 

25 ° 

86 ° 

0.18 

3.3 

4 . 0 ° 

1.07 

0 

2.29 

2 

35 

60 ° 

76 ° 

16 ° 

76 ° 

0.43 

2.2 

1 . 3 ° 

0.75 

1.20 

3 

33 

43 ° 

80 ° 

37 ° 

80 ° 

0.16 

2.5 

8 . 1 ° 

0.86 

1.50 

4 

33 

57 ° 

67 ° 

10 ° 

67 ° 

0.68 

1.4 

O 

CM 

• 

O 

0.30 

0.35 

10 

34 

22 ° 

49 ° 

39 o 

61 ° 

0.35 

1.4 

0 . 8 ° 

0.32 

0.38 

11 

34 

33 ° 

57 ° 

24 ° 

57 ° 

0.43 

1.5 

1 . 7 ° 

0.40 

0.50 

13 

34 

00 

o 

65 ° 

26 ° 

74 ° 

0.52 

1.9 

7 . 4 ° 

0.64 

0.94 

14 

34 

41 ° 

74 ° 

33 ° 

74 ° 

0.25 

2.8 

0 . 8 ° 

0.94 

1.76 

16 

35 

38 ° 

72 o 

33 ° 

71 ° 

0.26 

2.5 

1 . 2 ° 

0.87 

1.53 

18 

34 

55 ° 

77 * 

22 ° 

77 ° 

0.34 

2.1 

4 . 9 ° 

0.69 

1.06 

19 

34 

66 ° 

85 ° 

19 ° 

85 ° 

0.19 

3.6 

1 . 6 ° 

1.13 

2.58 

25 

35 

36 ° 

62 ° 

26 ° 

62 ° 

0.39 

1.7 

1 . 3 ° 

0.54 

0.74 

26 

35 

40 ° 

62 ° 

22 8 

62 ° 

0.45 

1.6 

O 

OO 

• 

O 

0.48 

0.63 

57 

41 

45 ° 

67 ° 

22 ° 

67 ° 

0.42 

1.8 

0 

OO 

• 

0 

0.58 

0.82 

62 

41 

70 ° 

82 ° 

12 8 

82 ° 

0.38 

2.5 

0 . 1 ° 

0.87 

1.53 

63 

41 

44 ° 

72 » 

27 ° 

71 ° 

0.32 

2.3 

0 

• 

O 

0.78 

1.29 

66 

41 

52 ° 

67 ° 

15 ° 

67 ° 

0.54 

1.6 

0.6° 

0.45 

0.59 

78 

41 

38 ° 

62 ° 

25 ° 

63 ° 

0.41 

1.7 

3 . 9 ° 

0.50 

0.67 
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and where a close choice had to be made among different S/C the event 
with the smallest 6 and/or largest a was chosen, where Judicious 
weighting between these rules was used on occasion. Only 18 shocks 
remain out of 22, as the table shows. Only two of the discarded cases 
(Codes 5 and 7) were due to the possibility of inaccurate normals 
(A = 23° and 68°, respectively) but most likely were due to poorly 
determined AB's. (Code 7 refers to a thick "shock" and undoubtedly had 
a rather poorly determined £,B but probably a good normal.) The other 2 
cases (Codes 80 and 82) arose because of small a. This was due to the 
shocks being nearly perpendicular types (8^ ~ 85°), and our interest in 
this section is in oblique shocks. Ideally P 2 = 8 ^ + Ot according to the 
magnetic coplanarity theorem. Comparing columns U and 6 shows that this 
closely holds in almost every case; one should note that this result was 
not forced (by a circular argument), since the associated normals for 
these shocks were almost exclusively obtained by multiple S/C methods. 

The table shows that in almost all cases R £ 1/2 and for this set 
a typical value is R = 1/3, the average being 0.37. The particular 
values of shown occur partly by selection, but it is interesting that 
nc rarely exceeds 35 s '— and recall that S/C data giving small a were 
discarded for this table. It is not obvious which of the three para- 
meters F./F, , S, or S best reflects probable shock strength, without 
2 1 o 

accompanying plasma data, but the obvious Inverse relationship between 
R and F^/Fj is intri 8 uin 8* Figure 8 shows a plot of Fj/F^ vs. R for 
these 18 shocks. 
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The figure suggests the statistical relationship (with a least 


squares rms - 0.39): 


F /F 

r i 


3.5 - 3.7 R 


where R ~ tan 0^/tan e ? . Chao (1970) demonstrates that 


F 2 sin B 2 


* 4- 1 


F x sin B l " M^OyNj)-! * 


(3) 


( 4 ) 


for a thermally isotropic medium, where is the pre-shock Allven Mach number 
and N is the plasma number density. Since the normal component of the 
magnetic field acrosc the shock front is continuous, l.e. since F 2 cos 


Fj cos 0 lt 


Eqs. 3 and 4 can je combined tc yield 

. f y F i + °- 2 

M A " J 3.7 » /N“ + F 2 /F 1 - 3.5 


(5) 


Hence, the strength of an oblique TP shock at 1 AU for the isotropic 
assumption depends, in a statistical way, only on the scalar quantities Fj/Fj 
and N^/N 2> In general the empirical relationship given by Eq.(5) may not be 
very accurate for any particular shock, but it suggests a typical property 
of oblique IP shocks a: 1 AU. However, we tested it for several previously 
studied specific shocks which were pa rametor-f itted according to schemes 
developed by Hhao (1970) and Lepplng and Argentiero (1971), and it yielded 
estimates of M. which when compared to best-fit values gave discrepancies 
ranging from 307. to only 107,. It appears that when N^/N^ gets too small 
(say < 0.3) relationship (5) yields unreliable estimates of M . 

Discussion 

Firstly, we would like to point out that our criterion for selecting 
SSC's, which was based on the agreement of reports from 20 or more magnetic 
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observatories, was consistent with assuming SSC events are caused by 
interplanetary shocks. From 48 events for which interplanetary magnetic 
field and/or plasma data were available, only 157. were caused by 
tangential discontinuities. For a continuation of the 9tudy described 
above and for completeness we relaxed our criterion for SSC-selection by 
including also those events identified by only a few (5 to 10) observatories. 

In those cases less than approximately 15% of the SSC's were caused by 

IP shocks; the Identification of these IP events was not always unambiguous. 

The results of the extended study are not reported here. Hence, large 

and universally observed SSC's tend overwhelmingly to be caused by IP 
shocks, and the lesser "SSC's" also may be caused by IP shocks but more 
likely by tangential discontinuities and other solar wind irhotnogeneltles. 

Secondly, we wish to stress the distinction between this statistical 
study and detailed studies on a single shock surface geometry (Mariani 
et al. 1970; Lepping and Chao, 1972). Figure 5 suggests that the average 
shock surface near the earth's orbit lies on a circle centered at the sun 
with a radius of 1 AU. However, the standard deviation from the mean 
surface is about 22°. This indicates that an individual shock surface 
geometry can deviate rather markedly from a spherical shape. From a 
physical point of view a spherically symmetric shoe’ near the sun propa- 
gating through IP space will Interact with interplanetary discontinuities 
(Lepping and Burlaga, 1973) and IP large scale inhomogeneities, such as 
streams and gradients in density, temperature, velocity, etc. (Heinemann 
and Slscoe, 1973). That is, it appears that an individual shock 
front may be distorted not only in gross geometry from a spherical shape. 
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but also experience a "rippling” on the scale of the correlation length 
of IP tangential discontinuities. A collection of these 
various perturbations on shock surfaces may result in an ensemble average 
of such surfaces being spherical with a relative large standard deviation 
even though few of the individual shock surfaces were actually spherical. 
In this sense the RMS deviation of the normals is probably more important 
than their average. 

Figures 5 and 6 (especially 5), showing a tendency for a radius of 
curvature of 1 AU at 1 AU regardless of the location of the flare site 
with respect to the shock observation point, tend to suggest that the 
initial (solar) shock shape is less important than IP processes in causing 
shock normal scattering over 1 AU (a % 23°). In fact, the symmetry 
indicated by these figures is rather remarkable from the viewpoint of IP 
shock normals. Itcwcver, an equally striking asymmetry exists which 
suggests a relationship about IP shock survivability. Notice that of 
the 33 shocks listed in Table 3 having a flare site association, 11 have 
an East solar longitude designation but 22 have a West designation. That 
is, it appears that an observer at 1 AU is twice as likely to find shocks 
related to West longitude flares as East longitude; the average for the 
set of 33 flares i* 17°W. Even the limited set of 22 shocks (Table 3, 
top, and Figure 5) give essentially the same results: 14 West longitude 

and 8 East longitude. Sakurai (1973 a,b) finds an apparently related 
asymmetrical tendency also based on SSC-flare associations: his 

statistical study shows that a maximum mean speed direction lies % 30%/ 
of the central meridian of the solar disk, i.e. near to the mean spiral 
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direction. The IP shock geometry for the Feb. 15-16, 1967 shock discussed 
by Lepping and Chao (1972) indicates a strong similarity. It appears 
that a shock experiences a greater chance of survival beyond 1 AU if it 
propagates along a direction approximately parallel to that average 
spiral which maps back to the flare site. This is not unreasonable 
according to the model by Heinemann and Siscoe (1973), which crucially 
depends on the IP spiral geometry and on large scale interactions of IP 
shocks with streams. The meridional plane statistics are less interesting: 
19 North latitude and 14 South latitude for the set of 33 shocks, and 13 
North and 9 South for the set of 22 shocks. 

In conclusion, interplanetary shock waves and geomagnetic storm 
sudden commencements (SSC's) are correlated with the solar activity of 
flares, radio bursts of Type II and/or Type IV. The average IP shock 
front at 1 AU has a radius of curvature on the order of 1 AU. However, 
the geometry of the shock front deduced to be of spherical shape is 
obviously not representative of any individual event, because an 
individual shock front may be severely distorted over 1 AU by IP processes. 
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Figure Captions 

Figure 1 Distribution of mean transit speeds over 1 AU computed using 

the difference in time between the occurrence of the flare and 


the onset of the SSC at earth for 81 events. 


Figure 2 The orientations of 38 shock normals. The angle is the solar 
ecliptic latitude and <f> s is the solar ecliptic longitude. The 
dotted arrows represent events with larger uncertainties than 
those with solid arrows (see text). 


Figure 3 W, the local shock speed in the solar wind frame versus the 

solar wind speed, n is an accurately estimated unit shock normal. 

Figure 4 A comparison of the local shock speed and the mean transit speed. 

Figure 5 Orientations in the ecliptic plane of the local shock surfaces at 
the appropriate heliocentric longitudes relative to their corres- 
ponding flares. The numbers associated with each event are the 
Code numbers (see Table 3). 


Figure 6 Orientations in the meridional plane of the local shock surfaces 
at the appropriate heliocentric latitudes relative to their 
corresponding flares. 

Figure 7 The distribution of mean transit speeds over 1 AU for simulated 
(random) cases (solid line) and for real cases (dashed line). 


The distribution of the real cases is the same as that in Figure 
1, but has been properly scaled so that the total number of 
events equals 46. 

Figure 8 F 2^ F 1 vs R tan 6^/tan 02» see Table 7) for 18 choice examples 

of flare associated shocks. 
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